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Abstract 

The tools of presymplectic geometry are used to study light rays trajectories in anisotropic 
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media. 
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The usefulness of the Lie-Hamilton optics in many different problems as ray tracing for 
ray design and computation of aberrations suggests the study of what happens for anisotropic 
media, because of the recent interest in the use of anisotropic optical material. This moti- 
vated a very recent paper |I| where the Hamiltonian formulation of geometric anisotropic 
Q_i| optics was studied. The theory was reexamined in ||. The basic principle of the theory is 

the celebrated Fermat's principle of least time (or extremal time if reflection is also allowed). 
In other words, light rays connecting points A and B are lines in the space in such a way 
that they satisfy the following variational condition: 

H ! rB 

5 n ds = 0. 



A 

The refractive index of the medium is given by the quotient n — -, and then physics tell 
us that n > 1. 

In a recent paper || we analysed from a geometric perspective the relationship for the 
case of isotropic media of the problem of determination of extremal curves for ([!]) with that 
of geodesies of a Riemannian metric conformal to the Euclidean metric. We will consider 
here the case in which the medium is not isotropic but it may depend on the velocity, or more 
specifically, on the direction of the ray. In this last case the problem cannot be reduced to a 
problem of Riemannian geometry as it happened when the refractive index n only depended 
on the position. So, the techniques of Presymplectic geometry are unavoidable for dealing 
with this dependence of the refractive index with the ray direction. 

To begin we remark the strong similarity of Fermat's principle with the more traditional 
Hamilton's principle of Classical Mechanics, with a Lagrangian function given by 



L = nyfg(v,v). (2) 
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This Lagrangian function is different iable only in the set of velocity phase space obtained 
by removing the null velocity points, i.e., the zero section of the tangent bundle. Moreover, 
the Lagrangian L is homogeneus of degree one, 

and consequently the corresponding energy function vanishes identically. Therefore the La- 
grangian is singular, because taking derivatives with respect to of both sides of the pre- 
ceding equation we obtain 

d 2 L 

dv l dvi ' ^ 

and then the Hessian matrix 

2 L 

is singular. The theory should be carefully reexamined using the tools of Presymplectic 
Geometry, as it was done in || for the isotropic media, where the tools of Riemannian 
Geometry were shown to be very useful. Actually, it is possible to show that in the latter case 
the solution curves for the regular Lagrangian system described by the regular Lagrangian 

^ = ^n 2 g(v,v) (6) 

are just the curves solution of the original problem, even if the curves are reparametrized 
(see e.g. [|J and references therein). Our aim is to analyse what happens in the more general 
case in which the refractive index can depend on the ray direction, i.e., the refractive index 
is a homogeneous function of degree zero of velocities in the set obtained by removing the 
zero section of the velocity phase espace. 

The geometric approach to Lagrangian Classical mechanics uses as velocity phase space 
the tangent bundle TM of the configuration space M that is assumed to be a differentiate 
manifold of dimension N. /,From now on we will follow the notation used in |J . The tangent 
structure is characterized by a (1, 1)— tensor field called vertical endomorphism S that in 
terms of natural coordinates (q l ,v l ) of the tangent bundle TM is given by 

S = ^ « (7) 

Given a function L e C°°(TM), we may define an exact 2-form in TM, ujl = —dOi, with 
the 1-form 9l being defined by 6l = dL o S, and a function El = A(L) — L, called energy 
function. In the above mentioned coordinates of TM we have the following expressions: 

a = < io > 
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Here A G X(TM) denotes the Liouville vector field generating dilations along the fibres. 
When lol is nondegenerate, i.e., the Hessian matrix (||) is regular, it defines a symplectic 
structure on TM, and a vector field Yl uniquely determined by i{Tl)^>l = dE^. 

We will next examine the problem of light rays, even for the more general case in which 
the refractive index n depends on the ray direction. This means that the refractive index 
must be homogeneous of degree zero in the velocities, An = 0. We define a new Lagrangian 
L = \L 2 . The Lagrangian L is homogeneous of degree two in the velocities, AIL = 2L = L 2 , 
and then E-n = L. 

Now taking into account that 6-^ = dh o S we see that the Liouville 1-form is 
proportional to 9l, namely, = L6l, (see ||) and then 

"L = jUJ h + -^dhA6 h . (12) 

As indicated above, when the refractive index n does not depend on velocities the 2-form 
tjjj^ is regular. In this more general case, however, it may be singular, because 

cj l = Luj l - dL A9 L (13) 

means that 

lu l AN = L N u L AN -N L N ~ l dLA9 L A co L A ^ = -N L N ^ dLA6 L A co L ^ N ~ 1 \ (14) 

and then u-^ AN can be identically null and in this case L would be singular. 

In the following we will restrict ourselves to the case in which L is regular and therefore 
there will be a uniquely defined vector field such that 

*(r L VL = dE h = d ^ (is) 

and then 

^(r L )(rfL) = r L (L) = o. (16) 

Moreover, is a second order differential equation vector field. 

First we check that the Liouville vector field A lies in keru;^. Indeed, 

z(AH = ^(AV L + Ah6 h -±dL(i(A)9 L ), (17) 

and taking into account that A is vertical and (h semibasic, the last term vanishes. Moreover, 

,(AV L = -C A 9^ + d(i(A)9jJ = -C A 9 h = -9 L (18) 
and therefore, taking into account that AIL = 2L = L 2 we find that 

i(A)u L = 0. (19) 
Secondly, is also in the kernel of u>l, because 
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and being a SODE, 



i(T L )# L = (dh o S)(T h ) = AL = 2L = L 2 



(21) 



and therefore 





and in particular 



t/( a Vl 



i 

1? 





(24) 



and therefore 




dim(kerco>£) = 2. 

Under these circumstances is possible to apply the reduction theory of presymplectic man- 
ifolds, following the ideas developed by Marsden and Weinstein ||. Presymplectic structures 
may arise either when using some constants of motion for reducing the phase space or also 
when the Lagrangian that has been chosen is singular. Then we will have a pair (Pq, Qq) 
where Qq is a closed but degenerate 2-form. A consistent solution of the dynamical equation 
can only be found in some points, leading in this way to the final constraint submanifold 
P introduced by Dirac (see e.g. ||). The pull back Q of the form Q on this manifold will 
be assumed to be of constant rank. The recipe for dealing with these systems was given by 
Marsden and Weinstein ||. First, in every point m G P, kerf2 m is a fc-dimensional linear 
space, so defining what is called a fc-dimensional distribution. The important point is that 
closedness of Q is enough to warrant that the distribution is integrable (and then it is called 
foliation): for any point m e P, there is a /c-dimensional submanifold of P passing through 
m and such that the tangent space at any point wl of this surface coincides with kerf2 m /. 
Such integral fc-dimensional submanifolds give a foliation of P by disjoint leaves and in the 
case in which the quotient space P = Pj ker Q is a differentiable manifold, then it is possible 
to define a nondegenerate closed 2-form Q in P such that tt*Q = Q. Here Jr : P — > P is the 
natural projection. It suffices to define Cl(vx,v 2 ) = fi(t>i,t>2), where V\ and v 2 are tangent 
vectors to P projecting under 7f* onto V\ and v 2 respectively. The symplectic space (P, Q) is 
said to be the reduced space. We will illustrate the method finding coordinates adapted to 
the distribution defined by the kernel ker u>l of the presymplectic structure defined by the 
singular optical Lagrangian in the case of a system in which either the index n depends on 
the third coordinate x 3 alone or the very interesting case in which the system is anisotropic 
and n is a function of the ray direction. We will determine the quotient reduced space and 
we will look for Darboux coordinates in this reduced symplectic manifold. Once Darboux 
coordinates have been found we can consider the problem from the active viewpoint and 
take advantage of the algebraic methods recently developed for computing aberrations (see 



e.g. 0). 
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Let us now consider the most general isotropic case in which the refractive index of the 
medium is not constant but it is given by a smooth function n(x 1 ,x 2 ,x 3 ). Fermat's princi- 
ple suggests us to consider the corresponding mechanical problem described by a singular 
Lagrangian L(q,v) = [g(v , v )] 1 ^ 2 , where g is a metric conformal to the Euclidean metric go? 

g(v,w) = n 2 g (v,w). (26) 



This problem was analysed in || where, as above is cited, it was shown that its study 
can be reduced to that of a regular Lagrangian L = \L 2 . This La grangian L is quadratic in 
velocities and the dynamical vector field solution of the dynamical equation i{TjJu^ = 
dE-^ = dh is not only a second order differential equation vector field but, moreover, it is a 



spray [JXTJ , the projection onto R 3 of its integral curves being the geodesies of the Levi-Civita 
connection defined by g. Then, is the geodesic spray given by 



r 



H =v i -^--T i jk v j v k ^- 
L dq l 3 dv* 



(27) 



where the Christoffel symbols T l aj. are 



jk 



dg k i , dgn dg 



+ 



-Jjk 



dqi dq k dq l 



(2* 



with g i: > being the inverse matrix of gij. 

In the particular case we are considering where g(v,w) = n 2 g (v,w), it was also shown 
above that the kernel of ujl is two-dimensional and it is generated by Tj^ and the Liouville 
vector field A. The distribution defined by kerc^L is integrable because ojl is closed; actually 
[A, TjJ = T]^ and the distribution is also generated by A and K defined by K = for 
which [A, K] = 0. In cartesian coordinates the Christoffel symbols are expressed as follows: 



jk 



1 

n 



dn ^ dn ^ dn ^ 
dxi k dx k j dx i k 



and the vector field K is given as 



K = ~ 

v 



—v l (v.Vn) 
n 



(29) 



dn \ d 



n dx l ) dv l 



(30) 



The theory of distributions suggests us the introduction of new local coordinates y l = 
F l (x, v), i = 1, . . . , 6, adapted to the distribution defined by ker cu L , i.e., such that K = 

A = (see |TI|). The search for these new coordinates is based on the solution of the 
partial differential equation system 

KF 1 = 1, AF 1 = 0, KF 2 = 0, AF 2 = 1, 

and 



KF 2+a _ g AF : 



2+a 



0, for a 



5 



The explicit computation of these functions depends very much on the choice of the 
function ^(x 1 , x 2 , x 3 ). We will illustrate next the theory with an particular example. If n 
only depends on x 3 , the presymplectic form can be written in the way 



dx 1 A d 



nv 



nv 



v l2 + v 22 + v 32 



+ 



n v 



32 



and the dynamical vector field is 



3/2 



+ dx 2 A d I = = (31) 

7 J 

J d ( ^ A rfx 3 + ^ ( ^ A ^3] ^ (32) 

^ <9x* n dx 3 dv 1 n dx 3 9w 2 ^ 

+ + ^ _^3 (34) 

n ax 13 air 

After some calculations we find the solution of the former systems (see |§). According 
to this, we will do the following choice for the new coordinates: 



y 



1 _ v _ 2 



x 2 - 



X 



C 3 



^(o-^Da+Ci 2 ) 



1 z 1 92 1 q2 
U + f + "IT 

2 



nv 



1 2 1 9 2 1 q2 



log 



n\l v l2 + f 22 + f 32 



where 



and 



c - 



C, = n 



\ 



v^ 2 + v* 2 



doing the inverse change and after some easy calculations remains 



ti L = d[y 1 + ^-y 2 ) Ady^ + dy 2 A dy 5 



which shows that 

e 



y 1 + -7V 2 = x l 



X 



v 2 J° V(" 2 (0-c 3 2 )(i + c?) 

d( 



y/(n 2 (C)-C$)(l + C?) 



(35) 

(36) 

(37) 
(38) 

(39) 
(40) 

(41) 
(42) 

(43) 

(44) 
(45) 
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and the corresponding 



1 4 

v = y 



nv 



v l2 + v 22 + v 32 ' 



v 2 = y 5 



nv 



^v l2 + v 22 + v* 2 



(46) 



are Darboux coordinates for the symplectic form induced in the quotient space. 

Let us now consider the particular but important case case in which the refractive index 
becomes constant out of a region. If for x 3 > L, the index n is constant, the above mentioned 
Darboux coordinates £ x and £ 2 are 



x — 



C 3 x 3 



rr 



x 2 - 



C 3 x 3 



y-cg)(i + c? 



up to a constant, and from the expresions of C\ and C 3 we see that the Darboux coordinates 
become 



1 1,1 3 



nv 



nv 



x -x 



v l2 + v 22 + v 32 



v l2 + v 22 + v 32 ' 



(47) 



in full agreement with ||13|| . Therefore, for an optical system such that the refractive index 
depends only on x 3 and, furthermore, the region in which the index is not constant is 
bounded, we can choose Darboux coordinates by fixing a a; 3 outside this region and taking 
Darboux coordinates for the corresponding problem of constant index. This justify the choice 
of coordinates as usually done for the ingoing and outgoing light rays in the corresponding 
constant index media, i.e. it shows the convenience of using flat screens in far enough regions 
on the left and right respectively, and then this change of Darboux coordinates seems to be, 



from an active viewpoint, a canonical transformation (see fT4||). 

We will next find the symplectic structure arising in an anisotropic optical medium, as 
well as some Darboux coordinates for it. We recall what we are only considering anisotropic 
media for which the refractive index depends only on the ray direction, i.e., n = n(v) but 
An = 0. In this case the presymplectic form remains as 



d 2 n Vi dn 
W\\^—n h 



+ 



dxWvi 
d 2 n 



Ml dxi 



dx l A dx J 



+ 



Vi dn Va dn 



dvjdvi \\v\\ dvj 



+ 



n 



n 



If II dvj 



dx % A dv 



31 



and the vector field associated with L 



We still have that 



r 



L 



A,r 



d_ 

' dx % 



JL 



(48) 
(49) 

(50) 
(51) 



and then the her ul defines an involutive, and hence an integrable, distribution that is also 
generated by A and K, K being the vector field 



k = — r 



L 



(52) 
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commuting with A. In this way, we can find new local coordinates adapted to the distribution 
that allow us to find later on the symplectic form in the quotient manifold, by solving the 
following differential equation systems: 



A/ 

A/ 
A/ 





1 





Kf 

Kf 
Kf 



0. 

0. 
1. 



(53) 

(54) 
(55) 



Acoording to the solution of the former systems we will do the following choice for the new 
coordinates: 



yi 



-z — X, 



V 



X 2 



— , V2 = — , 



-z - y, x 3 



y 3 = log v z 



the inverse change being given by 



V1X3 - x x , 
Vi exp y 3 , 



y 2 x 3 - x 2 , 
y 2 exp y 3 , 



x 3 

expy 3 . 



(56) 
(57) 



(58) 
(59) 



Moreover, in the reduction of the presymplectic form we must point out that the condition 
An = is written in the new coordinates as J^- = 0. Finally, using the change of coordinate 
former, we get after some calculations the following symplectic form in the quotient manifold 



d 



+ d 



nyi 



1 dft 

y\ + y\ + 1 OVI 



ny 2 



1 dft 

yj + yl + l 9y 2 



A dxi 
A dx 2 



(60) 
(61) 



which is in full agreement with 

As a final comment, let us remark that, even in this case, if we restrict ourselves to a 
region of constant index, we recover the Darboux coordinates for constant index medium 
and we can think of the relation betwen the ingoing and outgoing light rays as a change 
of Darboux coordinates. Therefore this change of Darboux coordinates seems to be again, 
from an active viewpoint, a canonical transformation. The transformations of phase space 
will be, in general, non lineal, i.e., they generate optical aberrations. It is possible, finally, 
to analyse these aberrations using both group theoretical and Lie algebraic tools (see []TTJ, 
[T3|j, and references therein). 
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